function y = func_chi_mean_from_omega(omega_vec,CONSTS,plot_data)

    num_iter=1e+5; %1000;
    chi_from_w_sum = zeros(size(omega_vec,1),1);
    chi_from_w = zeros(size(omega_vec,1),num_iter);
    for ni=1:num_iter

        chi_from_w(:,ni) = func_chi_from_omega(omega_vec,CONSTS,false);
        chi_from_w_sum(:) = chi_from_w_sum(:) + chi_from_w(:,ni);
    end
    
    y = chi_from_w_sum./num_iter;
    
    if(plot_data)
        figure; plot(omega_vec, y);
    end

end